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SYSTEM, METHOD, AND PROGRAM PRODUCT FOR RE- 
PARAMETERIZING THREE DIMENSIONAL MODELS REPRESENTED AS 
C ATMULL-CLARK SUBDIVISION SURFACES 

FIELD OF THE INVENTION 

The invention relates to the field of three dimensional computer modeling. More 
specifically, the invention relates to re-parameterization of three dimensional computer 
models. 

15 BACKGROUND OF THE INVENTION 

Subdivision surfaces are emerging as a powerful representation for shape design. Their 
simplicity and potential to overcome difficulties associated with traditional spline-based 
modeling has made them a popular choice for several applications, including industrial 
20 modeling, movie production, and games. 

Subdivision defines a smooth surface recursively, as the limit of a sequence of meshes. A 
finer mesh is obtained from a coarse mesh by applying a set of fixed refinement rules, for 
example the CatmuU-Clark rules (see CatmuU, E. and Clark, J., Recusively Generated B- 
25 Spline Surfaces on Arbitrary Topological Meshes, CAD 10(6), pp. 350-355, 1978). The 
initial, coarsest mesh is also known as the base mesh. Figure 2 shows a prior art example 
of a base mesh and a sequence of refinements that lead to a smooth surface. A 
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5 subdivision surface (Figure 2E) is generated from an initial control polyhedron (Figure 
2A) by repeated subdivision (Figures 2B-2D). Control meshes obtained at several 
intermediate subdivision steps are shown in the middle (Figures 2B-2D). In this example, 
the CatmuU-Clark subdivision scheme was applied. 

10 The initial coarsest mesh serves as a base domain over which the resulting surface is 
naturally parameterized. See Figure 3. Intuitively, the correspondence between control 
points on a given subdivision level and points in the base domain can be illustrated using 
two copies of the base mesh and applying in parallel midpoint subdivision on one copy 
(Figure 3 A) and CatmuU-Clark subdivision on the other (Figure 3B), as illustrated in 

15 Figure 3. This ability to cast subdivision surfaces as parametric surfaces over a domain is 
crucial for many operations on this type of representation. 

Some formal background on the definition and properties of natural parameterizations is 
now provided. A review of existing methods for exact evaluation at arbitrary parameter 
20 values is also included. 

Throughout this disclosure, we use the term control point (or simply point) to refer to 
points computed by subdivision. We also use the terai vertex to designate vertices of the 
polyhedron that serves as a base domain, as well as vertices added to it by midpoint 
25 subdivision. We also make use of the concepts of valence and extraordinary vertex / 
point. The valence of a vertex is an integer number equal to the number of faces incident 
to that vertex, A vertex is termed boundary vertex if it lies on the boundary of an open 
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5 mesh. A vertex which is not a boundary vertex is called an interior vertex. For Catmull- 
Clark meshes, an interior vertex of the base mesh which has a valence different than four 
or a boundary vertex with valence different than two is deemed extraordinary (see vertex 
marked with circle 210 in Figure 2A). The concepts of l-ring and 2-ring around a vertex 
allow us to refer to the mesh neighborhood aroxmd that vertex. The 1-ring around a vertex 

10 is defined as the set of vertices belonging to faces adjacent to that vertex. The 2-ring 
around a vertex is defined as the set formed by the union of its 1-ring and the 1 -rings of 
the vertices in the 1-ring different from the vertex itself. Parallel definitions can be 
obtained from the definitions above by replacing the word 'Vertex" with the word '"point" 
to define the concepts of boundary point, interior point, extraordinary point, and 1-ring 

15 and 2-ring around a point. 

More formally, let D denote the subset of corresponding to the base domain. 
Repeated subdivision generates a sequence of piecewise linear maps defined over the 
common domain, D . When this sequence converges uniformly (as it is the case for the 

20 Catmull-Clark scheme), the limit map / :D-> defines the limit surface of 

subdivision. For the Catmull-Clark scheme, the limit surface around an extraordinary 
point is completely defined by the set of points in the 2-ring around the extraordinary 
point. This set is called a control set. Let /[p] denote the limit surface generated by 
subdivision from the control setp . For converging subdivision schemes, the limit 

25 function of subdivision can be decomposed into a sum of basis functions: 

/[p](t/,v) = 2a,y;,(«,v) (1) 
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5 These so-called eigenbasis functions are the limit functions generated by the eigenvectors 
of the subdivision matrix S : /) := /[x J , where x, is the eigenvector of the subdivision 

matrix S corresponding to the eigenvalue A, : Sxj = A,X| . It has been shown in Zorin, D. 

Subdivision and Multiresolution Surface Representations. Pd.D. Thesis, Caltech, 1997 
that the eigenbasis functions satisfy the following scaling relation: 

10 /,(|,^) = V,(u,v) (2) 

It is possible to directly evaluate subdivision surfaces and their derivatives up to any 
order at arbitrary paramet^ values without the need for explicit subdivision, as proposed 
in Stam, J., Exact Evaluation of CatmuU-Clark Subdivision Surfaces at Arbitrary 
15 Parameter Values, Proceedings of Siggraph 98, pp. 395-404, 1998 and in Zorin, D. and 
Kristjansson, D., Evaluation of Piecewise Smooth Subdivision Surfaces, Visual 
Computer, 2002. 

The evaluation algorithm described in Zorin, D. and Kristjansson, D., Evaluation of 
20 Piecewise Smooth Subdivision Surfaces, Visual Computer, 2002 is now reviewed in the 
context of the Catmull-Clark subdivision scheme. 

The evaluation problem can be stated as follows: "Given a quadrilateral face J of the 
base mesh parameterized over the unit square [0,1] x [0,1] and parametric coordinates 
25 (w, v) within this face, find the value of the subdivision surface / (m, v, y ) e at this 
parameter value." 
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If the given quadrilateral face contains no extraordinary vertices, the surface can be 
evaluated as a bi-cubic spline patch defined by the control vertices of the face after an 
additional subdivision step. If the face contains extraordinary vertices, we can assxime 
without restricting the generality of the problem, that it has a single extraordinary vertex 
10 (otherwise the problem can be reduced to this case after one subdivision step). 

Let y = (Uy V, j) be a point in the domain of the surface different from an extraordinary 
vertex. An important observation is that there exists a sufficiently fine subdivision level 
m such that y is located in a quadrilateral Q on level m and the control mesh of Q does not 
15 contain any extraordinary vertices. This means that the surface over g is polynomial and 
can be evaluated explicitly (including at pointy ) using bi-cubic spline evaluation. 

Obviously, for points closer to the extraordinary vertex, the level m is higher than for 
points farther away. This interpretation leads to a conceptual decomposition of the 
neighborhood around an extraordinary vertex into concentric layers. These layers 

20 correspond to the number of subdivision st^s required to isolate the location being 
evaluated from the extraordinary vertex. Formally, if 5 denotes the subdivision matrix, 
the control meshes of quadrilaterals on level m can be computed from the initial control 
mesh by a multiplication with the/w -th power of the subdivision matrix: S"" . These 
observations lead to a direct evaluation algorithm as outlined in the pseudo-code of 

25 Algorithm 1 and as illustrated in the flowchart of Figure 5. 
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PROBLEMS WITH THE PRIOR ART 



By definition, a surface / :D-^R^is -continuous if for any point x € D there is a 

neighborhood of jc , qD and a regular parameterization of / (U^ ) over the unit disk in 

the plane. Intuitively, this means that a surface is smooth aroxmd a given point if it 
10 behaves like a plane in its immediate vicinity. 

According to this definition, the proof that the limit surface generated by a subdivision 
scheme is C^-continuous around an extraordinary point amoimts to showing that a regular 
parameterization as required in the definition exists at that point. It has been shown that, 

15 in the case of CatmuU-Clark surfaces, such a parameterization exists and hence, they are 
C^-continuous (see below also). However, it turns out that at the extraordinary points, the 
natural parameterization discussed above, do not satisfy the regularity requirement for 
valences greater than 4. To prove that this is the case, consider the neighborhood around 
an extraordinary vertex and expressions (1) and (2). For simplicity, we assume that the 

20 extraordinary vertex is situated at the origin (i.e., a,. = 0 ). 

If we derivate expression (1) with respect tOM , we obtain: 

|^/[p](«,v) = Xa,|^/;.(i/,v) 

The expression of the derivative with respect to v is similar and omitted for brevity. 
25 From (2), we have: 
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5 and by induction, 

^fii—.—) = (2X.r^Mu,v) (3) 
dw ' 2"" 2'" du ' ' ^ ^ ^ 

From (3), it follows that for m cx> the value — (0,0) is unbounded if (2A,. ) > 1 . For the 

du 

CatmuU-Clark scheme, the subdominant eigenvalue is larger than Vi for valences greater 
than 4, which proves our claim that derivatives computed based on the natural 
10 parameterization diverge around extraordinary points corresponding to extraordinary 
vertices of valence greater than 4. Using the same argument, if (2>l. ) < 1 the value 

of — (0,0) is zero for m ^ oo . This is the case for extraordinary points corresponding to 
du 

extraordinary vertices of valence 3. Figure 3 illustrates the behavior of the iso-parameter 
lines aroxmd extraordinary vertices of valence 3 (Figure 3 A), 5 (Figure 3B), and 13 
15 (Figure 3C), respectively. Intuitively, the behavior of the derivatives can be observed by 
examining the spacing between iso-parameter lines. For valences greater than 4, the 
spacing between consecutive lines increases as they approach the extraordinary point as 
the derivatives diverge. For valence 3 the spacing between iso-parameter lines decreases 
as the derivatives converge to zero. 



20 



All references cited above are herein incorporated by reference in their entirety. 



25 
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5 ASPECTS OF THE INVENTION 



An aspect of the invention is an improved system and method for general three- 
dimensional surface modeling. 

10 An aspect of the invention is an improved system and method for re-parameterizing 
tiiree-dimensional models represented as Catmull-Clark subdivision surfaces. 

An aspect of the invention is an improved system and method for surface evaluation in 
soUd modeling. 

15 

SUMMARY OF THE INVENTION 

The present invention is a system and method for re-parameterizing one or more three- 
20 dimensional Catmull-Clark models. Each of the models has one or more surfaces and 
one or more extraordinary vertices. One or more adjacent iso-parameter lines of the 
model have a natural spacing between them that changes as the iso-parameter lines 
approach the extraordinary vertex. The invention includes a re-parameterization process 
that re-parameterizes the model so that one or more of the extraordinary vertices have 
25 adjacent iso-parameter lines with new spacing that is different than as the lines approach 
the extraordinary vertex. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a block diagram of one preferred system comprising the present invention. 

10 Figure 2, comprising Figures 2A-2E, illustrates a prior art subdivision process. 

Figure 3, comprising Figures 3 A and 3B, illustrates a prior art natural parameterization of 
subdivision surfaces. 

15 Figure 4, comprising Figures 4A-4C, illustrates the prior art behavior of the iso-parameter 
lines around extraordinary vertices of valence 3 (Figure 4 A), 5 (Figure 4B), and 13 
(Figure 4C), respectively. 

Figure 5 is a flow chart of a prior art surface evaluation process. 

20 

Figure 6 is a flow chart of a novel re-parameterization process. 

Figure 7 is a flow chart of one preferred vanishing derivative re-parameterization process 
that is used in the re-parameterization process of Figure 6. 

25 

Figure 8 is a flow chart of one preferred characteristic map re-parameterization process 
that is used as an alternative process in the re-parametmzation process of Figure 6. 
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Figure 9, comprising Figures 9A and 9B, illustrates the behavior of the iso-parameter 
lines around extraordinary vertices of valence 5 (Figure 9 A), and 13 (Figure 9B), 
respectively after re-parameterization with vanishing derivatives. 

10 Figure 10, comprising Figure 1 OA- IOC, illustrates prior examples of control meshes and 
corresponding characteristic maps for the Catmull-Clark scheme for valences 3 (Figure 
lOA), 4 (Figure lOB), and 5 (Figure IOC). 

Figure 1 1 shows an example prior art structure of layers and patches of the characteristic 
15 map around an extraordinary vertex. 

Figure 12, comprising Figures 12A-12C, shows iso-parameter lines corresponding to 
parameterization by inversion of tiie characteristic map around vertices of valence 3 
(Figure 12A), 5 (Figure 12B), and 13 (Figure 12C), respectively. 

20 

Figure 13, comprising Figures 13A and 13B, shows a graph and shape of a weighting 
function w for two different values of the parameter^ : ^ = 3 (Figure 1 3 A) and k = 6 
(Figure 13B). 

25 
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5 DETAIL DESCRIPTION OF THE INVENTION 



One aspect this disclosure addresses is the issue of robust, constant-time evaluation of 
subdivision surfaces and their first-order derivatives near extraordinary vertices. It is 
well-known that the evalxiation of first-order derivatives based on the natural 

10 parameterization of subdivision sxirfaces is unstable as the derivatives diverge near 
extraordinary vertices of valence greater than four. Two alternative novel preferred 
embodiments of re-parameterization methods around extraordinary vertices are 
presented. The first one ensures that the derivatives after re-parameterization converge to 
zero around extraordinary vertices. This method is relatively simple to implement, 

15 however vanishing derivatives may be problematic in certain applications. The second 
approach, although slightly more involved, provides a general robust method to compute 
a regular parameterization around extraordinary vertices with guaranteed well-behaved 
derivatives. It is based on the inversion of the characteristic map. 



20 Figure 1 is one example and simplified block diagram of a preferred data processing 
system 100 that is suitable for practicing this invention. The data processing system 100 
includes at least one data processor 101 coupled to a bus 102 through which the data 
processor 101 may address a memory sub-system 103, also referred to herein simply as 
the memory 103. The memory 103 may include RAM, ROM and fixed and removable 

25 disks and/or tape. The memory 103 is assumed to store a program containing program 
instructions for causing the data processor 101 to execute methods in accordance with the 
teachings of this invention. Also stored in the memory 103 can be at least one database 
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5 104 of digital 3D models, such as models of objects, a non-limiting example of which is 
shown in Figure 2. The digital model data may include models obtained by scanning real 
objects using a 3D scanner, digital models created using a computer, models imported 
from existing databases, and/or models generated by or derived from any suitable means. 
In general, the digital model data may be any desired type or types of 3D models, 

10 including but not limited to digital models of persons, places, animals, plants, mechanical 
parts and assemblies, packaging, containers, abstract forms and so forth. For the purpose 
of this invention, the models are assumed to be represented as Catmull-Clark subdivision 
surfaces with or without multiresolution details. 

15 The data processor 101 is also coupled through the bus 102 to a user interface, preferably 
a graphical user interface (GUI) 105 that includes any one or more of the following well 
known user interface devices: a keyboard, a mouse, a trackball, a voice recognition 
interface, and any general purpose user display device, such as a high resolution 
graphical CRT display terminal, a LCD display terminal, or any suitable other display 

20 device. 

The data processor 101 may also be coupled through the bus 102 to a network interface 
106 that provides bidirectional access to a data communications network 107, such as an 
intranet and/or the intemet. Coupled to the network 107 can be one or more sources 
25 and/or repositories of digital models, such as a remote digital model database 108 that is 
reachable through an associated server 109. 
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5 The data processor 101 is also preferably coupled through the bus 102 to at least one 
peripheral device 110, such as a scanner 1 lOA (e.g., a 3D scanner) and/or a printer HOB 
and/or a 3D model making apparatus, such as a rapid prototyping system, and/or a 
computer controlled fabrication system. 

10 In general, this invention may be implemented using one or more software programs 
running on a personal computer, a server, a microcomputer, a mainframe computer, a 
portable computer, and embedded computer, or by any suitable type of programmable 
data processor 101. 

15 Given this disclosure, many other system architectures would become known to one 
skilled in the art. 

The use of this invention substantially improves the evaluation of 3D model data 
represented as a CatmuU-Clark subdivision surface for many applications that require 

20 surface evaluation (for example, surface editing and boolean operations). The methods 
may be used to process the digital 3D model data stored in the 3D model database 104 
and/or in the remotely stored 3D model database 108 over the network 107 and in 
cooperation with the server 109. As but one example, a 3D input subdivision surface 
model to be evaluated could be remotely stored in the 3D model database 108a, while the 

25 evaluated limit sinface corresponding to the model could be stored in the local 3D model 
database 104. 
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5 The invention includes an evaluation process 600, shown in more detail in Figure 6. 

Figure 5 is a flow chart of a prior art surface evaluation process. 

A typical prior art algorithm for evaluating subdivision surfaces at arbitrary parameter 
10 values is as follows: 

Algorithm 1 P = Eval(Surface^u,vj) 

Given a location y = (w,v,y) in the domain of the Surface, step 510 comprises 
determining a layer number m of (w, v) , i.e., a sufficiently fine subdivision level m such 
1 5 that y is located in a quadrilateral Q on level m , and the control mesh of this quadrilateral 
does not contain any extraordinary vertices. Inside g the surface is polynomial and it can 
be evaluated explicitly at (m, v) xising bi-cubic spline evaluation (520) to obtain the 
evaluated point P on the surface. 

20 While these are typical main steps in prior art methods, their specifics may differ. For 
two different examples see Stam, J., Exact Evaluation of Catmull-Clark Subdivision 
Surfaces at Arbitrary Parameter Values, Proceedings of Siggraph 98, pp. 395-404, 1998 
and Zorin, D. and Kristjansson, D., Evaluation of Piecewise Smooth Subdivision 
Surfaces, Visual Computer, 2002. 

25 

Figure 6 is a flow chart of a novel re-parameterization process. 
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5 Since the limit surface in the case of the Catmull-Clark scheme is C^-continuous, we 
know that a regular parameterization around each extraordinary vertex exists. In order to 
implement an evaluator that generates derivatives that converge at extraordinary vertices, 
we have to find such re-parameterizations around these vertices. The general format of a 
new evaluation algorithm we propose involves using a re-parameterization mapping ;r as 

10 follows: 

Algorithm 2 P- EvalRparam(Surface, u,vj) 

In order to determine a new location for a given parametric point that has undergone a 
transformation due to a re-parameterization, we compute the image of the parametric 
15 point being evaluated through the re-parameterization mapping;r : (w',v',y) = {n:{u,v)J) 
(610) and use the newly computed parametric point v',7) as input to the prior art 
evaluation process 500 (also described in Algorithm 1). The resulting point on the surface 
is the output of the latter: P = Eval(Surface,u,vj) 

20 We propose two novel methods for generating a re-parameterization mapping 7t . 

The first one (see Figure 7) entails a simple mapping that ensures that the surface has 
continuous first order derivatives that converge to 0 around extraordinary points. Such a 
mapping is relatively easy to implement, however vanishing derivatives may cause 
25 problems for certain applications. The second approach (see Figure 8) is to use the 
inverse of the characteristic map for re-parameterization. Although it requires a more 
involved implementation, this method offers the most robust alternative backed by 
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theoretical guarantees that ensure a regular parameterization and hence, well-behaved 
first-order derivatives. We discuss these two methods next. 



The first re-parameterization method is described next. Figure 7 is a flow chart of one 
preferred vanishing derivative re-parameterization process that is used in the re- 
10 parameterization process of Figure 6. 



F 1 

Our approach is to use a re-parameterization of the form: x > | x |^ " x around 

each vertex of a base mesh face j , where we denote by x = (m, v) the parametric point 
being evaluated inside face j . At each vertex, according to (3), in order to ensure that the 
15 derivatives converge after re-parameterization we need to impose the condition 

log A 

where /I is the subdominant eigenvalue of the subdivision matrix corresponding to the 
valence of that vertex. 



20 The idea behind the choice of constant in the right-hand side of expression (4) is to 

induce a scaling relation of the form (3), but with a scaling factor that converges to zero 
as m tends to infinity. If such a factor is chosen to be, for instance, of the 
form2*''^ A witha as in (4), then convergence to zero is ensured. 

25 Thus, we first must compute the four subdominant eigenvalues of the subdivision 

matrices corresponding to the four vertices of face j in which evaluation takes place 
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5 (710). Then we select appropriate vahxesaj^.k = 1,...,4 which satisfy condition (4) (720) 

ensuring that, after a re-parameterization of the form x > I x | * x , the 

derivatives aroxmd vertex A: converge to zero, for all A: = 1,,..,4 .If any of the vertices have 
valence less than or equal to four we use the natural parameterization around those 
vertices, i.e., oCf^ =1. 

10 

Using this kind of re-parameterization around the vertices of the base domain we ensure 
that the function defining the surface has continuous derivatives that converge to zero 
around those vertices. 

15 We can now create a re-parameterization mapping;ro (730) by blending together the four 
re-parameterization functions ,k = 1,...,4 using a weighted simi as defined in (5); 

Xw(|r,(i/,v)|)r,(i^, (?;(«, V))) 

7r(u,v) = -^ (5) 

J,H\T,(u,v)\) 

where (x) =| x |^*~^ x is the re-parameterization at vertex k, T^is a rotation fiinction 

20 that places the extraordinary vertex at the origin of the parameter space, and wis a 
weighting fiinction satisfying boundary conditions: 

wiO) = 1, w(l) = 0, w(r) = w\r) = 0, r = 0,1 . 

For examples of weighting functions see section 3.3. 
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Figure 9, comprising Figures 9A-9B, illustrates the behavior of the iso-parameter lines 
around extraordinary vertices of valence 5 (Figure 9 A), and 13 (Figure 9B), respectively 
after re-parameterization with vanishing derivatives. Note that as the valence increases, 
10 the spacing between consecutive lines decreases. 

The second method for generating a re-parameterization mapping;? makes use of the 
concept of a characteristic map. This notion is briefly reviewed next. 

15 It has been shown that control points near a vertex converge, up to an affine 

transformation, to a limit configuration that is independent of the control points of the 
original mesh (see Zorin, D. Subdivision and Multiresolution Surface Representations. 
Pd.D. Thesis, Caltech, 1997). This result motivates the introduction of the concept of a 
characteristic map. Informally, any subdivision surface generated by a given subdivision 

20 scheme behaves near an extraordinary vertex of valence k like the characteristic map of 
that scheme for valence k. Formally the characteristic map is defined as follows: 

Definition The characteristic map of a subdivision scheme for valence k is the map 
O : f/ ^ defined as O = fiol^i^x + ^2^2] ^ where e^ and e^ are unit coordinate 
25 vectors, Xy and jCj are subdominant eigenvectors, and /2/>[p] : t/ is the limit 
surface of subdivision generated by subdivision on U from the 2D control points p . 



YOR920030614US1 



18 



5 Figure 10, comprising Figure lOA-lOC, illustrates prior examples of characteristic maps 
for the CatmuU-Clark scheme for valences 3 (Figure lOA), 4 (Figure lOB), and 5 (Figure 
IOC). The control mesh in each of these figures is shown as 101 0 and images of 
regularly sampled points through the characteristic map are shown as 1000. For valence 
5, the structure of layers and patches around the extraordinary vertex is also illustrated in 
10 Figure IOC as 1020. 

The characteristic map plays an important role in the analysis of smoothness properties of 
subdivision schemes. 

15 The CatmuU-Clark scheme along with other popular subdivision schemes has been 

proven to be & -continuous at extraordinary vertices. For such schemes, the characteristic 
map O can be inverted and a parameterization of the surface around the extraordinary 
vertex of the form: / (O"' (s)) , where s e 0(f/) can be obtained. This observation stands 
at the basis of our new re-parameterization scheme. 

20 

Figure 8 is a flow chart of one preferred characteristic map re-parameterization process 
that is used as an altemative process in the re-parameterization process of Figure 6. 

F 

We use a re-parameterization of the form: x ><E>"^ (x) around each vertex of a base 

25 mesh face j , where we denote by x = (m, v) the point being evaluated inside face J . 
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5 First we compute the four characteristic mapsOj ,...,€>4 corresponding to the four vertices 
of face j in which evaluation takes place (810). Then we invert the characteristic maps 
obtaining <I)"S,.,.,0"\ (820) and we use these inverses to re-parameterize the region 
around each of the four vertices using a re-parameterization of the 
F 

formx 

10 

Using this kind of re-parameterization around the vertices of the base domain we ensure 
that the function defining the surface has continuous derivatives that converge to the 
actual derivative of the surface aroimd those vertices. 

1 5 We can now create a re-parameterization mapping Tt^^ (830) by blending together the four 
re-parameterization functions using a weighted simi as defined in (5), with 
(x) = (x) being the re-parameterization at vertex k. 

The characteristic map defines a smooth variation of parametric positions in the 
20 parameter domain around the extraordinary vertex. By re-parameterizing the surface 
around the extraordinary vertex using this map, we ensure that a similarly smooth 
variation occurs on the surface. 

We avoid the need to explicitly compute an inverse of the characteristic map by taking 
25 advantage of an organization into layers and patches around the extraordinary vertex as 
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5 shown in Figures IOC and Figure 1 1 . The numbers in Figure 1 1 indicate an example 
indexing scheme of patches within a layer. 

To compute the image of a point (m, v) in the domain of the surface through the inverse 
map we use the following procedure: 

10 

Algorithm 3 («>') = O"^ (w, v) 

First we determine the layer L of the characteristic map and the polynomial patch 
within L that contains the point (w,v) . Layer! corresponds to a sufficiratly fine 
subdivision level such that(M, v) is located in a patch on that. Inside we can use 
15 inverse polynomial evaluation to obtain (w\v*) . 

As shown in Figure 1 1, the domain is bounded by parametric curves rather than straight 
lines. To locate the layer and patch containing a given point in this type of domain we 
first determine an approximate location of the given point (t/,v) using bounding boxes of 

20 layers and patches. We then refine the search for the exact layer and patch by searching 
in a restricted region around the approximate location. For this we need to classify the 
position of the point (w,v) with respect to the curves bounding the layers and patches, as 
shown in Figure 1 1 . We transform these cubic polynomial curves firom their parametric 
representation to an implicit one using a process of imphcitization. Implicitization refers 

25 to the process of extracting an implicit representation for a curve given in parametric 
form. Such a problem arises often when attempting to decide whether a given point lies 
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5 on a given curve or not. The implicit representation offers an efficient way of finding the 
answer through a simple comparison test. 

One example method for implicitzation is described in Sederberg, T.W. and Chen, F., 
Implicitization Using Moving Curves and Surfaces, Proceedings of SIGGRAPH'95, pp. 
10 301-308, 1995. An alternative method is based on computing resultants and is briefly 
described next. 

Definition. Let p(t) = a^x"' +a^_^x""^ +... + aoand^(0 = ^„-x" +i>„_i^""^ +... + Z?o denote 
two scalar polynomials such that at least one of and6„ is non-zero. The Sylvester 

IS resultant of the two polynomials is defined as the determinant of the following 

(m + «) X (m + n) Sylvester matrix: 



«0 ^1 ^2 

0 ao 



0 0 



0 




0 0 0 



... 6, ... b„ 



20 The link between the Sylvester resultant and the common roots of the polynomials is 
given by the following well-known result: 
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5 Theorem. Polynomials p and q have a common root if and only if det S{p, q)-0. 



We use the Sylvester resultants to solve the implicitization problem. 

The implicitization of a parametric polynomial curve C{t) = (x(0, y(t)) using 
10 polynomial resultants is described next. 

Let (Xq , € Cbe an arbitrary point on the curve and let us consider the 

polynomials p^(t) = x(t) - and p^ (t) - y(t) - . Clearly, since (jCq, j^^) is on the 

curve, a common root tQ of the polynomials p^ and Py must exist. According to the 

previous theorem, this can happen if and only if det S(p^^Py) = 0 . Thus, 

15 dQtS(p^,p ) = 0 constitutes the implicit representation of the curve C . 



20 



In our case, the curves delineating layers and patches are defined by cubic polynomials 
(see Figure 11): 

Hence, the implicit representation of such curves is extremely simple: 



r(x,j;):det 



0 
0 

A -y 

0 
0 



0 



Br 



D, 0 

c, A 



5. 



A,-y B 



0 



2 '-2 

A2-y B, 



c, 

0 



0 
0 

A 

0 
0 



= 0 
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5 To determine on what side of the curve r(jc, y) a given point ( jCq , y^^ ) is, we simply 
compare r(xQ , j;q ) to 0 (or some £ to account for numerical precision errors). 

Once a patch containing the input point (m, v) is found, we must compute the pre-image of 
this point through the inverse of the patch. While there are several options to compute the 
10 inverse of a polynomial surface patch, we illustrate here one such method based on 
Newton's approach to root finding. However, we note that any polynomial patch 
inversion method could be used instead. 



Given a polynomial patch(P(uW) and a point (m, v) inside the patch, we need to find a 
15 point (w',vO such liiat(uW) = ^'^ {u,v) . In other words, we have to find (i/', v^) that 
satisfies the equation: 



= 0 (6) 



Newton's method is a root finding algorithm that uses the first few terms of the Taylor 
series of a fimction in the vicinity of a suspected root to approximate the actual. 

20 

If Tlu.v") = (fP, {u\v\T^ Wy)) is the patch containing (m,v) , let 
Q{ u\v) = (iPj (w', vO - w,lp2 («', v') - v)be the cubic polynomial corresponding to equation 
(6). We are looking for a root of the equation g(tt',v') = 0 inside the domain [0,1] x [0,1] . 
To find it, we apply Newton-type iterations, starting firom an initial guess at the center of 
25 the domain: 
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K,v;) = (o.5A5) 

We iterate until the approximation error becomes smaller than a given threshold. In our 
experiments we have observed the convergence to be very fast, typically requiring less 
than five iterations. 



10 Figure 12, comprising Figures 12A-12C, shows iso-parameter lines corresponding to re- 
parameterization by inversion of the characteristic map aroimd vertices of valence 3 
(Figure 12A), 5 (Figure 12B), and 13 (Figure 12C), respectively. Note that as the valence 
increases, the spacing between consecutive lines remains constant. 



15 Each of the two re-parameterization methods described in this invention makes use of a 
blending procedure to accommodate base faces with more than one extraordinary vertex. 
To blend re-parameterizations aroxmd different vertices of the same base face we use 
weighting functions that smoothly fall off away firom the vertex. An example choice of 
weighting function has the form 

20 

w(r) = B(kr-l)'hBikr) + B{kr + l) (7) 
with k a positive integer constant and ^ the piecewise cubic B-spline function defined as: 
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B(t) = 



- + 2t+t^ +-t\-2<t^-l 
3 6 

3 2 

3 2 

--2? + ^^ --^^l<^<2 
3 6 
0,2<t 



(8) 



The graphs of this function for values of ^ = 3 and A: = 6 are shown in Figure 13. 
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